home *** CD-ROM | disk | FTP | other *** search
/ Amiga Format CD 46 / Amiga Format CD46 (1999-10-20)(Future Publishing)(GB)[!][issue 1999-12].iso / -in_the_mag- / reader_requests / scilab / tests / sparse.dia.ref < prev    next >
Text File  |  1999-09-16  |  36KB  |  772 lines
  1.  
  2.  
  3. mode(5)
  4. lines(0)
  5. ij=[1 6;1 5;1 3;2 4;2 1;4 4;4 3;5 1;6 6];
  6. v=[1;2;3;4;5;6;7;8;9];
  7. ij1=[1 6;1 2;6 5];
  8. v1=[-1;-1;-1];
  9. vc=v+%i*(21:29)';
  10. v1c=v1+%i*[0;0.3;-1.2];
  11. zer=sparse([],[],[6,6]);
  12. //
  13. a=sparse([1 6;1 5;1 3;2 4;2 1;4 4;4 3;5 1;6 6],1:9,[6 6]);
  14. b=sparse([1 6;1 2;6 5],[-1;-1;-1],[6 6]);
  15. ma=sparse([1 6;1 5;1 3;2 4;2 1;4 4;4 3;5 1;6 6],-(1:9),[6 6]);
  16. mb=sparse([1 6;1 2;6 5],-[-1;-1;-1],[6 6]);
  17. ac=sparse([1 6;1 5;1 3;2 4;2 1;4 4;4 3;5 1;6 6],(1:9)+%i*(21:29),[6 6]);
  18. bc=sparse([1 6;1 2;6 5],[-1;-1;-1]+%i*[0;0.3;-1.2],[6 6]);
  19. mac=sparse([1 6;1 5;1 3;2 4;2 1;4 4;4 3;5 1;6 6],-(1:9)-%i*(21:29),[6 6]);
  20. mbc=sparse([1 6;1 2;6 5],[1;1;1]+%i*[0;-0.3;1.2],[6 6]);
  21. //
  22. if %t then
  23. //-------------------------------------------------------------
  24. //test des primitives sparse, full
  25. //--------------------------------------------------------------
  26. //  -- for real matrix
  27. if or(full(sparse(0.3))<>0.3) then bugmes();quit;end
  28. v=0*ones(1,3);v(3)=5;v(1)=1;
  29. if or(full(sparse([1 1;1 3],[1 5]))<>v) then bugmes();quit;end
  30. v=0*ones(1,300);v(300)=5;v(1)=1;
  31. if or(full(sparse([1 1;1 300],[1 5]))<>v) then bugmes();quit;end
  32. v=0*ones(3,1);v(3)=5;v(1)=1;
  33. if or(full(sparse([1 1;3 1],[1 5]))<>v) then bugmes();quit;end
  34. v=0*ones(300,1);v(300)=5;v(1)=1;
  35. if or(full(sparse([1 1;300 1],[1 5]))<>v) then bugmes();quit;end
  36. sp=sparse([1 6;1 5;1 3;2 4;2 1;4 4;4 3;5 1;6 6],1:9,[6 6]);
  37. a=[0,0,3,0,2,1;;5,0,0,4,0,0;0,0,0,0,0,0;0,0,7,6,0,0;8,0,0,0,0,0;0,0,0,0,0,9];
  38. if or(full(sp)<>a) then bugmes();quit;end
  39. sp=sparse([1 6;1 5;1 3;2 4;2 1;4 4;4 3;5 1;6 6],1:9,[8 6]);a(8,6)=0;
  40. if or(full(sp)<>a) then bugmes();quit;end
  41. if or(full(sparse([],[],[4 10]))<>0*ones(4,10)) then bugmes();quit;end
  42. v=sparse([2 1;3 1;4 1;6 1],[10;11;12;13],[6,1]);
  43. a=[0;10;11;12;0;13];
  44. if or(full(v)<>a) then bugmes();quit;end
  45. v=sparse([1 2;1 3;1 4;1 6],[10;11;12;13],[1,6]);
  46. a=[0,10,11,12,0,13];
  47. if or(full(v)<>a) then bugmes();quit;end
  48. //  -- for complex matrix
  49. sp=sparse([1 6;1 5;1 3;2 4;2 1;4 4;4 3;5 1;6 6],(1:9)+%i*(21:29),[6 6]);
  50. a=[0,0,3+%i*23,0,2+%i*22,1+%i*21
  51.    5+%i*25,0,0,4+%i*24,0,0
  52.    0,0,0,0,0,0
  53.    0,0,7+%i*27,6+%i*26,0,0
  54.    8+%i*28,0,0,0,0,0
  55.    0,0,0,0,0,9+%i*29];
  56. if or(full(sp)<>a) then bugmes();quit;end
  57. sp=sparse([1 6;1 5;1 3;2 4;2 1;4 4;4 3;5 1;6 6],(1:9)+%i*(21:29),[8 6]);
  58. a(8,6)=0;
  59. if or(full(sp)<>a) then bugmes();quit;end
  60. v=sparse([2 1;3 1;4 1;6 1],[10-3*%i;11;12+5*%i;13+0.5*%i],[6,1]);
  61. a=[0;10-%i*3;11;12+%i*5;0;13+%i*0.5];
  62. if or(full(v)<>a) then bugmes();quit;end
  63. v=sparse([1 2;1 3;1 4;1 6],[10-3*%i;11;12+5*%i;13+0.5*%i],[1,6]);
  64. a=[0,10-%i*3,11,12+%i*5,0,13+%i*0.5];
  65. if or(full(v)<>a) then bugmes();quit;end
  66. //-------------------------------------------------------------
  67. //test de la primitives spget
  68. //--------------------------------------------------------------
  69. sp=sparse([1 6;1 5;1 3;2 4;2 1;4 4;4 3;5 1;6 6],1:9,[6 6]);
  70. [ij,v]=spget(sp);
  71. if or(ij<>[1,3;1,5;1,6;2,1;2,4;4,3;4,4;5,1;6,6])|or(v<>[3;2;1;5;4;7;6;8;9]) then bugmes();quit;end
  72. zer=sparse([],[],[6,6]);[ij,v]=spget(zer);
  73. if ij<>[]|v<>[] then bugmes();quit;end
  74. //-----------------------------------------------
  75. // addition and substraction tests
  76. //-----------------------------------------------
  77. // -- real real
  78. a=sparse([1 6;1 5;1 3;2 4;2 1;4 4;4 3;5 1;6 6],1:9,[6 6]);
  79. b=sparse([1 6;1 2;6 5],[-1;-1;-1],[6 6]);
  80. v=sparse([2 1;3 1;4 1;6 1],[10;11;12;13],[6,1]);
  81. vt=sparse([1 2;1 3;1 4;1 6],[10;11;12;13],[1,6]);
  82. ma=sparse([1 6;1 5;1 3;2 4;2 1;4 4;4 3;5 1;6 6],-(1:9),[6 6]);
  83. mb=sparse([1 6;1 2;6 5],-[-1;-1;-1],[6 6]);
  84. zer=sparse([],[],[6,6]);
  85. ac=sparse([1 6;1 5;1 3;2 4;2 1;4 4;4 3;5 1;6 6],(1:9)+%i*(21:29),[6 6]);
  86. bc=sparse([1 6;1 2;6 5],[-1;-1;-1]+%i*[0;0.3;-1.2],[6 6]);
  87. vc=sparse([2 1;3 1;4 1;6 1],[10-3*%i;11;12+5*%i;13+0.5*%i],[6,1]);
  88. vct=sparse([1 2;1 3;1 4;1 6],[10-3*%i;11;12+5*%i;13+0.5*%i],[1,6]);
  89. mac=sparse([1 6;1 5;1 3;2 4;2 1;4 4;4 3;5 1;6 6],-(1:9)-%i*(21:29),[6 6]);
  90. mbc=sparse([1 6;1 2;6 5],[1;1;1]+%i*[0;-0.3;1.2],[6 6]);
  91. if or(full(v+sparse([],[],[6,1]))<>full(v)) then bugmes();quit;end
  92. if or(full(sparse([],[],[6,1])+v)<>full(v)) then bugmes();quit;end
  93. if or(full(v+v)<>full(v)+full(v)) then bugmes();quit;end
  94. if or(full(v-v)<>full(v)-full(v)) then bugmes();quit;end
  95. if or(full(vt+sparse([],[],[1,6]))<>full(vt)) then bugmes();quit;end
  96. if or(full(vt+vt)<>full(vt)+full(vt)) then bugmes();quit;end
  97. if or(full(vt-vt)<>full(vt)-full(vt)) then bugmes();quit;end
  98. if or(full(zer+zer)<>0*ones(6,6)) then bugmes();quit;end
  99. if or(full(a+a)<>full(a)+full(a)) then bugmes();quit;end
  100. if or(full(b+b)<>full(b)+full(b)) then bugmes();quit;end
  101. if or(full(a+zer)<>full(a)) then bugmes();quit;end
  102. if or(full(zer+a)<>full(a)) then bugmes();quit;end
  103. if or(full(b+a)<>full(b)+full(a)) then bugmes();quit;end
  104. if or(full(a+b)<>full(b)+full(a)) then bugmes();quit;end
  105. if or(full(a+ma)<>full(a)+full(ma)) then bugmes();quit;end
  106. if or(full(a-a)<>full(a)-full(a)) then bugmes();quit;end
  107. if or(full(a-ma)<>full(a)-full(ma)) then bugmes();quit;end
  108. if or(full(b-mb)<>full(b)-full(mb)) then bugmes();quit;end
  109. if or(full(a-zer)<>full(a)) then bugmes();quit;end
  110. if or(full(zer-a)<>-full(a)) then bugmes();quit;end
  111. if or(full(a-mb)<>full(a)-full(mb)) then bugmes();quit;end
  112. //  -- real complex
  113. if or(full(sparse([],[],[6,1])+vc)<>full(vc)) then bugmes();quit;end
  114. if or(full(v+vc)<>full(v)+full(vc)) then bugmes();quit;end
  115. if or(full(v-vc)<>full(v)-full(vc)) then bugmes();quit;end
  116. if or(full(vt+vct)<>full(vt)+full(vct)) then bugmes();quit;end
  117. if or(full(vt-vct)<>full(vt)-full(vct)) then bugmes();quit;end
  118. if or(full(a+ac)<>full(a)+full(ac)) then bugmes();quit;end
  119. if or(full(b+bc)<>full(b)+full(bc)) then bugmes();quit;end
  120. if or(full(a+bc)<>full(a)+full(bc)) then bugmes();quit;end
  121. if or(full(zer+ac)<>full(zer)+full(ac)) then bugmes();quit;end
  122. if or(full(b+ac)<>full(b)+full(ac)) then bugmes();quit;end
  123. if or(full(a-ac)<>full(a)-full(ac)) then bugmes();quit;end
  124. if or(full(b-bc)<>full(b)-full(bc)) then bugmes();quit;end
  125. if or(full(a-bc)<>full(a)-full(bc)) then bugmes();quit;end
  126. if or(full(zer-ac)<>full(zer)-full(ac)) then bugmes();quit;end
  127. if or(full(b-ac)<>full(b)-full(ac)) then bugmes();quit;end
  128. // -- complex real
  129. if or(full(vc+v)<>full(vc)+full(v)) then bugmes();quit;end
  130. if or(full(vc-v)<>full(vc)-full(v)) then bugmes();quit;end
  131. if or(full(vct+vt)<>full(vct)+full(vt)) then bugmes();quit;end
  132. if or(full(vct-vt)<>full(vct)-full(vt)) then bugmes();quit;end
  133. if or(full(ac+a)<>full(ac)+full(a)) then bugmes();quit;end
  134. if or(full(bc+b)<>full(bc)+full(b)) then bugmes();quit;end
  135. if or(full(ac+b)<>full(ac)+full(b)) then bugmes();quit;end
  136. if or(full(ac+zer)<>full(zer)+full(ac)) then bugmes();quit;end
  137. if or(full(bc+a)<>full(bc)+full(a)) then bugmes();quit;end
  138. if or(full(ac-a)<>full(ac)-full(a)) then bugmes();quit;end
  139. if or(full(bc-b)<>full(bc)-full(b)) then bugmes();quit;end
  140. if or(full(ac-b)<>full(ac)-full(b)) then bugmes();quit;end
  141. if or(full(ac-zer)<>-full(zer)+full(ac)) then bugmes();quit;end
  142. if or(full(bc-a)<>full(bc)-full(a)) then bugmes();quit;end
  143. // -- complex complex
  144. if or(full(vc+vc)<>full(vc)+full(vc)) then bugmes();quit;end
  145. if or(full(vc-vc)<>full(vc)-full(vc)) then bugmes();quit;end
  146. if or(full(vct+vct)<>full(vct)+full(vct)) then bugmes();quit;end
  147. if or(full(vct-vct)<>full(vct)-full(vct)) then bugmes();quit;end
  148. if or(full(ac+ac)<>full(ac)+full(ac)) then bugmes();quit;end
  149. if or(full(bc+bc)<>full(bc)+full(bc)) then bugmes();quit;end
  150. if or(full(ac+bc)<>full(ac)+full(bc)) then bugmes();quit;end
  151. if or(full(bc+ac)<>full(bc)+full(ac)) then bugmes();quit;end
  152. if or(full(ac-ac)<>full(zer)) then bugmes();quit;end
  153. if or(full(bc-bc)<>full(zer)) then bugmes();quit;end
  154. if or(full(ac-bc)<>full(ac)-full(bc)) then bugmes();quit;end
  155. if or(full(bc-ac)<>full(bc)-full(ac)) then bugmes();quit;end
  156. // sparse full and full sparse operation (soft coded apoerations)
  157. if or(full(a+1)<>full(a)+1) then bugmes();quit;end
  158. if or(full(1+a)<>full(a)+1) then bugmes();quit;end
  159. if or(full(a+ones(6,6))<>full(a)+ones(6,6)) then bugmes();quit;end
  160. if or(full(ones(6,6)+a)<>full(a)+ones(6,6)) then bugmes();quit;end
  161. if or(full(a+2*eye)<>full(a)+2*eye) then bugmes();quit;end
  162. if or(full(2*eye+a)<>full(a)+2*eye) then bugmes();quit;end
  163. if or(full(a-1)<>full(a)-1) then bugmes();quit;end
  164. if or(full(1-a)<>1-full(a)) then bugmes();quit;end
  165. if or(full(a-ones(6,6))<>full(a)-ones(6,6)) then bugmes();quit;end
  166. if or(full(ones(6,6)-a)<>ones(6,6)-full(a)) then bugmes();quit;end
  167. if or(full(a-2*eye)<>full(a)-2*eye) then bugmes();quit;end
  168. if or(full(2*eye-a)<>2*eye-full(a)) then bugmes();quit;end
  169. if or(full(ac+1)<>full(ac)+1) then bugmes();quit;end
  170. if or(full(1+ac)<>full(ac)+1) then bugmes();quit;end
  171. if or(full(ac+ones(6,6))<>full(ac)+ones(6,6)) then bugmes();quit;end
  172. if or(full(ones(6,6)+ac)<>full(ac)+ones(6,6)) then bugmes();quit;end
  173. if or(full(ac+2*eye)<>full(ac)+2*eye) then bugmes();quit;end
  174. if or(full(2*eye+ac)<>full(ac)+2*eye) then bugmes();quit;end
  175. if or(full(ac-1)<>full(ac)-1) then bugmes();quit;end
  176. if or(full(1-ac)<>1-full(ac)) then bugmes();quit;end
  177. if or(full(ac-ones(6,6))<>full(ac)-ones(6,6)) then bugmes();quit;end
  178. if or(full(ones(6,6)-ac)<>ones(6,6)-full(ac)) then bugmes();quit;end
  179. if or(full(ac-2*eye)<>full(ac)-2*eye) then bugmes();quit;end
  180. if or(full(2*eye-ac)<>2*eye-full(ac)) then bugmes();quit;end
  181. if or(full(ac+full(bc))<>full(ac)+full(bc)) then bugmes();quit;end
  182. if or(full(ac-full(bc))<>full(ac)-full(bc)) then bugmes();quit;end
  183. if or(full(full(ac)+full(bc))<>full(ac)+full(bc)) then bugmes();quit;end
  184. end
  185. if %t then
  186. //-----------------------------------------------
  187. // multiplication tests
  188. //-----------------------------------------------
  189. // real x real
  190. // sparse scalar , saclar sparse
  191. if or(full(a*2)<>full(a)*2) then bugmes();quit;end
  192. if or(full(2*a)<>full(a)*2) then bugmes();quit;end
  193. if  a*[]<>[] then bugmes();quit;end
  194. if  []*a<>[] then bugmes();quit;end
  195. c=rand(6,6);
  196. if norm(a*c-full(a)*c) >100*%eps then bugmes();quit;end
  197. if norm(c*a-c*full(a)) >100*%eps then bugmes();quit;end
  198. // sparse sparse
  199. if or(full(zer*zer)<>full(zer)) then bugmes();quit;end
  200. if or(full(a*zer)<>full(zer)) then bugmes();quit;end
  201. if or(full(zer*a)<>full(zer)) then bugmes();quit;end
  202. if  norm(full(a*a)-full(a)*full(a),1)>100*%eps then bugmes();quit;end
  203. if  norm(full(b*b)-full(b)*full(b),1)>100*%eps then bugmes();quit;end
  204. if  norm(full(a*b)-full(a)*full(b),1)>100*%eps then bugmes();quit;end
  205. if  norm(full(b*a)-full(b)*full(a),1)>100*%eps then bugmes();quit;end
  206. // complex x real real x complex
  207. // sparse scalar , scalar sparse
  208. if or(full(ac*2)<>full(ac)*2) then bugmes();quit;end
  209. if or(full(2*ac)<>full(ac)*2) then bugmes();quit;end
  210. if norm(full((2+%i)*a)-(2+%i)*full(a),1) >100*%eps then bugmes();quit;end
  211. if norm(full(a*(2+%i))-(2+%i)*full(a),1) >100*%eps then bugmes();quit;end
  212. if  ac*[]<>[] then bugmes();quit;end
  213. if  []*ac<>[] then bugmes();quit;end
  214. c=rand(6,6);
  215. cc=c+rand(6,6)*%i;
  216. if norm(ac*c-full(ac)*c) >100*%eps then bugmes();quit;end
  217. if norm(cc*a-cc*full(a)) >100*%eps then bugmes();quit;end
  218. // sparse sparse
  219. if or(full(ac*zer)<>full(zer)) then bugmes();quit;end
  220. if or(full(zer*ac)<>full(zer)) then bugmes();quit;end
  221. if  norm(full(ac*a)-full(ac)*full(a),1)>100*%eps then bugmes();quit;end
  222. if  norm(full(bc*b)-full(bc)*full(b),1)>100*%eps then bugmes();quit;end
  223. if  norm(full(a*bc)-full(a)*full(bc),1)>100*%eps then bugmes();quit;end
  224. if  norm(full(b*ac)-full(b)*full(ac),1)>100*%eps then bugmes();quit;end
  225. // // complex x complex
  226. if norm(ac*cc-full(ac)*cc) >100*%eps then bugmes();quit;end
  227. if norm(cc*ac-cc*full(ac)) >100*%eps then bugmes();quit;end
  228. // sparse sparse
  229. if  norm(full(ac*ac)-full(ac)*full(ac),1)>100*%eps then bugmes();quit;end
  230. if  norm(full(bc*bc)-full(bc)*full(bc),1)>100*%eps then bugmes();quit;end
  231. if  norm(full(bc*ac)-full(bc)*full(ac),1)>100*%eps then bugmes();quit;end
  232. //----------------------------------------------------------
  233. // element wise multiplication tests
  234. //----------------------------------------------------------
  235. if or(full(ac.*2)<>full(ac)*2) then bugmes();quit;end
  236. if or(full((2).*ac)<>full(ac)*2) then bugmes();quit;end
  237. if  a.*[]<>[] then bugmes();quit;end
  238. if  [].*a<>[] then bugmes();quit;end
  239. c=rand(6,6);
  240. //if norm(a*c-full(a)*c) >100*%eps then bugmes();quit;end
  241. //if norm(c*a-c*full(a)) >100*%eps then bugmes();quit;end
  242. // sparse sparse
  243. if or(full(zer.*zer)<>full(zer)) then bugmes();quit;end
  244. if or(full(a.*zer)<>full(zer)) then bugmes();quit;end
  245. if or(full(zer.*a)<>full(zer)) then bugmes();quit;end
  246. if  norm(full(a.*a)-full(a).*full(a),1)>100*%eps then bugmes();quit;end
  247. if  norm(full(b.*b)-full(b).*full(b),1)>100*%eps then bugmes();quit;end
  248. if  norm(full(a.*b)-full(a).*full(b),1)>100*%eps then bugmes();quit;end
  249. if  norm(full(b.*a)-full(b).*full(a),1)>100*%eps then bugmes();quit;end
  250. // complex x real real x complex
  251. // sparse scalar , scalar sparse
  252. if norm(full((2+%i).*a)-(2+%i).*full(a),1) >100*%eps then bugmes();quit;end
  253. if norm(full(a.*(2+%i))-(2+%i).*full(a),1) >100*%eps then bugmes();quit;end
  254. if  ac.*[]<>[] then bugmes();quit;end
  255. if  [].*ac<>[] then bugmes();quit;end
  256. c=rand(6,6);
  257. cc=c+rand(6,6)*%i;
  258. if norm(full(ac.*c)-full(ac).*c) >100*%eps then bugmes();quit;end
  259. if norm(full(cc.*a)-cc.*full(a)) >100*%eps then bugmes();quit;end
  260. // sparse sparse
  261. if or(full(ac.*zer)<>full(zer)) then bugmes();quit;end
  262. if or(full(zer.*ac)<>full(zer)) then bugmes();quit;end
  263. if  norm(full(ac.*a)-full(ac).*full(a),1)>100*%eps then bugmes();quit;end
  264. if  norm(full(bc.*b)-full(bc).*full(b),1)>100*%eps then bugmes();quit;end
  265. if  norm(full(a.*bc)-full(a).*full(bc),1)>100*%eps then bugmes();quit;end
  266. if  norm(full(b.*ac)-full(b).*full(ac),1)>100*%eps then bugmes();quit;end
  267. // // complex x complex
  268. if norm(ac.*cc-full(ac).*cc) >100*%eps then bugmes();quit;end
  269. if norm(cc.*ac-cc.*full(ac)) >100*%eps then bugmes();quit;end
  270. // sparse sparse
  271. if  norm(full(ac.*ac)-full(ac).*full(ac),1)>100*%eps then bugmes();quit;end
  272. if  norm(full(bc.*bc)-full(bc).*full(bc),1)>100*%eps then bugmes();quit;end
  273. if  norm(full(bc.*ac)-full(bc).*full(ac),1)>100*%eps then bugmes();quit;end
  274. // ----------------------------------------------------------
  275. // test de la transposition
  276. //-----------------------------------------------------------
  277. if or(full(a')<>full(a)') then bugmes();quit;end
  278. if or(full(ac')<>full(ac)') then bugmes();quit;end
  279. if or(full(zer')<>full(zer)' ) then bugmes();quit;end
  280. v=sparse([2 1;3 1;4 1;6 1],[10;11;12;13],[6,1]);
  281. if or(full(v')<>full(v)') then bugmes();quit;end
  282. if or(full((v')')<>full(v)) then bugmes();quit;end
  283. vc=sparse([2 1;3 1;4 1;6 1],[10-3*%i;11;12+5*%i;13+0.5*%i],[6,1]);
  284. if or(full(vc')<>full(vc)') then bugmes();quit;end
  285. if or(full((vc')')<>full(vc)) then bugmes();quit;end
  286. // ----------------------------------------------------------
  287. // test des concatenation
  288. //-----------------------------------------------------------
  289. if or(full([a])<>full(a)) then bugmes();quit;end
  290. if or(full([a b])<>[full(a) full(b)]) then bugmes();quit;end
  291. if or(full([a;b])<>[full(a);full(b)]) then bugmes();quit;end
  292. if or(full([a []])<>full(a)) then bugmes();quit;end
  293. if or(full([a;[]])<>full(a)) then bugmes();quit;end
  294. if or(full([a zer])<>[full(a) full(zer)]) then bugmes();quit;end
  295. if or(full([zer;b])<>[full(zer);full(b)]) then bugmes();quit;end
  296. if or(full([ac])<>full(ac)) then bugmes();quit;end
  297. if or(full([ac b])<>[full(ac) full(b)]) then bugmes();quit;end
  298. if or(full([ac;b])<>[full(ac);full(b)]) then bugmes();quit;end
  299. if or(full([ac []])<>full(ac)) then bugmes();quit;end
  300. if or(full([ac;[]])<>full(ac)) then bugmes();quit;end
  301. if or(full([a bc])<>[full(a) full(bc)]) then bugmes();quit;end
  302. if or(full([a;bc])<>[full(a);full(bc)]) then bugmes();quit;end
  303. if or(full([ac bc])<>[full(ac) full(bc)]) then bugmes();quit;end
  304. if or(full([ac;bc])<>[full(ac);full(bc)]) then bugmes();quit;end
  305. // ----------------------------------------------------------
  306. // test des extractions
  307. //-----------------------------------------------------------
  308. af=full(a);
  309. if or(full(a(1,3))<>af(1,3)) then bugmes();quit;end
  310. if or(full(a(1,4))<>af(1,4)) then bugmes();quit;end
  311. if or(full(a(1,:))<>af(1,:)) then bugmes();quit;end
  312. v=sparse([2 1;3 1;4 1;6 1],[10;11;12;13],[6,1]);vf=full(v);
  313. if or(full(v(:))<>vf(:)) then bugmes();quit;end
  314. if or(full(v(3:4))<>vf(3:4)) then bugmes();quit;end
  315. if or(full(v([1 5]))<>vf([1 5])) then bugmes();quit;end
  316. if or(full(v([4 3]))<>vf([4 3])) then bugmes();quit;end
  317. if or(full(v([4 4]))<>vf([4 4])) then bugmes();quit;end
  318. if or(full(v([1 1]))<>vf([1 1])) then bugmes();quit;end
  319. v=v';vf=vf';
  320. if or(full(v(:))<>vf(:)) then bugmes();quit;end
  321. if or(full(v(3:4))<>vf(3:4)) then bugmes();quit;end
  322. if or(full(v([1 5]))<>vf([1 5])) then bugmes();quit;end
  323. if or(full(v([4 3]))<>vf([4 3])) then bugmes();quit;end
  324. if or(full(v([4 4]))<>vf([4 4])) then bugmes();quit;end
  325. if or(full(v([1 1]))<>vf([1 1])) then bugmes();quit;end
  326. acff=full(ac);
  327. if or(full(ac(1,3))<>acff(1,3)) then bugmes();quit;end
  328. if or(full(ac(1,4))<>acff(1,4)) then bugmes();quit;end
  329. if or(full(ac(1,:))<>acff(1,:)) then bugmes();quit;end
  330. vc=sparse([2 1;3 1;4 1;6 1],[10-3*%i;11;12+5*%i;13+0.5*%i],[6,1]);vcf=full(vc);
  331. if or(full(vc(:))<>vcf(:)) then bugmes();quit;end
  332. if or(full(vc(3:4))<>vcf(3:4)) then bugmes();quit;end
  333. if or(full(vc([1 5]))<>vcf([1 5])) then bugmes();quit;end
  334. if or(full(vc([4 3]))<>vcf([4 3])) then bugmes();quit;end
  335. if or(full(vc([4 4]))<>vcf([4 4])) then bugmes();quit;end
  336. if or(full(vc([1 1]))<>vcf([1 1])) then bugmes();quit;end
  337. vc=vc';vcf=vcf';
  338. if or(full(vc(:))<>vcf(:)) then bugmes();quit;end
  339. if or(full(vc(3:4))<>vcf(3:4)) then bugmes();quit;end
  340. if or(full(vc([1 5]))<>vcf([1 5])) then bugmes();quit;end
  341. if or(full(vc([4 3]))<>vcf([4 3])) then bugmes();quit;end
  342. if or(full(vc([4 4]))<>vcf([4 4])) then bugmes();quit;end
  343. if or(full(vc([1 1]))<>vcf([1 1])) then bugmes();quit;end
  344. // ----------------------------------------------------------
  345. // test des insertions
  346. //-----------------------------------------------------------
  347. end
  348. a=sparse([1 6;1 5;1 3;2 4;2 1;4 4;4 3;5 1;6 6],1:9,[6 6]);
  349. vt=sparse([1 2;1 3;1 4;1 6],[10;11;12;13],[1,6]);
  350. // full line insertion
  351. //----------------------
  352. a1=a;a1(1,:)=vt;A=full(a);Vt=full(vt);A(1,:)=Vt;
  353. if or(full(a1)<>A) then bugmes();quit;end
  354. a1=a;a1(3,:)=vt;A=full(a);Vt=full(vt);A(3,:)=Vt;
  355. if or(full(a1)<>A) then bugmes();quit;end
  356. a1=a;a1(1,:)=vt;A=full(a);Vt=full(vt);A(1,:)=Vt;
  357. if or(full(a1)<>A) then bugmes();quit;end
  358. a1=a;a1(3,:)=vt;A=full(a);Vt=full(vt);A(3,:)=Vt;
  359. if or(full(a1)<>A) then bugmes();quit;end
  360. a1=a;a1(7,:)=vt;A=full(a);Vt=full(vt);A(7,:)=Vt;
  361. if or(full(a1)<>A) then bugmes();quit;end
  362. a1=a;a1(8,:)=vt;A=full(a);Vt=full(vt);A(8,:)=Vt;
  363. if or(full(a1)<>A) then bugmes();quit;end
  364. a1=a;a1([1 3],:)=[vt;2*vt];A=full(a);Vt=full(vt);A([1 3],:)=[Vt;2*Vt];
  365. if or(full(a1)<>A) then bugmes();quit;end
  366. a1=a;a1([3 1],:)=[vt;2*vt];A=full(a);Vt=full(vt);A([3 1],:)=[Vt;2*Vt];
  367. if or(full(a1)<>A) then bugmes();quit;end
  368. a1=a;a1([1 3 1],:)=[vt;2*vt;3*vt];A=full(a);Vt=full(vt);A([1 3 1],:)=[Vt;2*Vt;3*Vt];
  369. if or(full(a1)<>A) then bugmes();quit;end
  370. //  insert zero vector
  371. vt=sparse([],[],[1,6]);
  372. a1=a;a1(1,:)=vt;A=full(a);Vt=full(vt);A(1,:)=Vt;
  373. if or(full(a1)<>A) then bugmes();quit;end
  374. a1=a;a1(3,:)=vt;A=full(a);Vt=full(vt);A(3,:)=Vt;
  375. if or(full(a1)<>A) then bugmes();quit;end
  376. a1=a;a1(1,:)=vt;A=full(a);Vt=full(vt);A(1,:)=Vt;
  377. if or(full(a1)<>A) then bugmes();quit;end
  378. a1=a;a1(3,:)=vt;A=full(a);Vt=full(vt);A(3,:)=Vt;
  379. if or(full(a1)<>A) then bugmes();quit;end
  380. a1=a;a1(7,:)=vt;A=full(a);Vt=full(vt);A(7,:)=Vt;
  381. if or(full(a1)<>A) then bugmes();quit;end
  382. a1=a;a1(8,:)=vt;A=full(a);Vt=full(vt);A(8,:)=Vt;
  383. if or(full(a1)<>A) then bugmes();quit;end
  384. a1=a;a1([1 3],:)=[vt;2*vt];A=full(a);Vt=full(vt);A([1 3],:)=[Vt;2*Vt];
  385. if or(full(a1)<>A) then bugmes();quit;end
  386. a1=a;a1([3 1],:)=[vt;2*vt];A=full(a);Vt=full(vt);A([3 1],:)=[Vt;2*Vt];
  387. if or(full(a1)<>A) then bugmes();quit;end
  388. a1=a;a1([1 3 1],:)=[vt;2*vt;3*vt];A=full(a);Vt=full(vt);A([1 3 1],:)=[Vt;2*Vt;3*Vt];
  389. if or(full(a1)<>A) then bugmes();quit;end
  390. a=sparse([],[],[6,6]);
  391. a1=a;a1(1,:)=vt;A=full(a);Vt=full(vt);A(1,:)=Vt;
  392. if or(full(a1)<>A) then bugmes();quit;end
  393. a1=a;a1(3,:)=vt;A=full(a);Vt=full(vt);A(3,:)=Vt;
  394. if or(full(a1)<>A) then bugmes();quit;end
  395. a1=a;a1(1,:)=vt;A=full(a);Vt=full(vt);A(1,:)=Vt;
  396. if or(full(a1)<>A) then bugmes();quit;end
  397. a1=a;a1(3,:)=vt;A=full(a);Vt=full(vt);A(3,:)=Vt;
  398. if or(full(a1)<>A) then bugmes();quit;end
  399. a1=a;a1(7,:)=vt;A=full(a);Vt=full(vt);A(7,:)=Vt;
  400. if or(full(a1)<>A) then bugmes();quit;end
  401. a1=a;a1(8,:)=vt;A=full(a);Vt=full(vt);A(8,:)=Vt;
  402. if or(full(a1)<>A) then bugmes();quit;end
  403. b=sparse([1 1;1 3;1 6;2 1;2 2;2 4],10:15,[2,6]);
  404. a1=a;a1([1 3],:)=b;A=full(a);B=full(b);A([1 3],:)=B;
  405. if or(full(a1)<>A) then bugmes();quit;end
  406. a1=a;a1([3 1],:)=b;A=full(a);B=full(b);A([3 1],:)=B;
  407. if or(full(a1)<>A) then bugmes();quit;end
  408. b=sparse([1 1;1 3;1 6;2 1;2 2;2 4;3 3;3 5],10:17,[3,6]);
  409. a1=a;a1([1 3 1],:)=b;A=full(a);B=full(b);A([1 3 1],:)=B;
  410. if or(full(a1)<>A) then bugmes();quit;end
  411. //  insert zero vector
  412. vt=sparse([],[],[1,6]);
  413. a1=a;a1(1,:)=vt;A=full(a);Vt=full(vt);A(1,:)=Vt;
  414. if or(full(a1)<>A) then bugmes();quit;end
  415. a1=a;a1(3,:)=vt;A=full(a);Vt=full(vt);A(3,:)=Vt;
  416. if or(full(a1)<>A) then bugmes();quit;end
  417. a1=a;a1(1,:)=vt;A=full(a);Vt=full(vt);A(1,:)=Vt;
  418. if or(full(a1)<>A) then bugmes();quit;end
  419. a1=a;a1(3,:)=vt;A=full(a);Vt=full(vt);A(3,:)=Vt;
  420. if or(full(a1)<>A) then bugmes();quit;end
  421. a1=a;a1(7,:)=vt;A=full(a);Vt=full(vt);A(7,:)=Vt;
  422. if or(full(a1)<>A) then bugmes();quit;end
  423. a1=a;a1(8,:)=vt;A=full(a);Vt=full(vt);A(8,:)=Vt;
  424. if or(full(a1)<>A) then bugmes();quit;end
  425. b=sparse([1 1;1 3;1 6;2 1;2 2;2 4],10:15,[2,6]);
  426. a1=a;a1([1 3],:)=b;A=full(a);B=full(b);A([1 3],:)=B;
  427. if or(full(a1)<>A) then bugmes();quit;end
  428. a1=a;a1([3 1],:)=b;A=full(a);B=full(b);A([3 1],:)=B;
  429. if or(full(a1)<>A) then bugmes();quit;end
  430. b=sparse([1 1;1 3;1 6;2 1;2 2;2 4;3 3;3 5],10:17,[3,6]);
  431. a1=a;a1([1 3 1],:)=b;A=full(a);B=full(b);A([1 3 1],:)=B;
  432. if or(full(a1)<>A) then bugmes();quit;end
  433. // full column insertion
  434. //----------------------
  435. a=sparse([1 6;1 5;1 3;2 4;2 1;4 4;4 3;5 1;6 6],1:9,[6 6]);
  436. v=sparse([2 1;3 1;4 1;6 1],[10;11;12;13],[6,1]);
  437. a1=a;a1(:,1)=v;A=full(a);V=full(v);A(:,1)=V;
  438. if or(full(a1)<>A) then bugmes();quit;end
  439. a1=a;a1(:,2)=v;A=full(a);V=full(v);A(:,2)=V;
  440. if or(full(a1)<>A) then bugmes();quit;end
  441. a1=a;a1(:,1)=v;A=full(a);V=full(v);A(:,1)=V;
  442. if or(full(a1)<>A) then bugmes();quit;end
  443. a1=a;a1(:,2)=v;A=full(a);V=full(v);A(:,2)=V;
  444. if or(full(a1)<>A) then bugmes();quit;end
  445. a1=a;a1(:,3)=v;A=full(a);V=full(v);A(:,3)=V;
  446. if or(full(a1)<>A) then bugmes();quit;end
  447. //
  448. a1=a;a1(:,7)=v;A=full(a);V=full(v);A(:,7)=V;
  449. if or(full(a1)<>A) then bugmes();quit;end
  450. a1=a;a1(:,8)=v;A=full(a);V=full(v);A(:,8)=V;
  451. if or(full(a1)<>A) then bugmes();quit;end
  452. b=sparse([1 2;2 1;3 1;3 2;4 1;6 2],10:15,[6,2]);
  453. a1=a;a1(:,[1 3])=b;A=full(a);B=full(b);A(:,[1 3])=B;
  454. if or(full(a1)<>A) then bugmes();quit;end
  455. a1=a;a1(:,[3 1])=b;A=full(a);B=full(b);A(:,[3 1])=B;
  456. if or(full(a1)<>A) then bugmes();quit;end
  457. b=sparse([1 2;2 1;2 3;3 1;3 2;4 1;5 3;6 2],10:17,[6,3]);
  458. a1=a;a1(:,[1 3 1])=b;A=full(a);B=full(b);A(:,[1 3 1])=B;
  459. if or(full(a1)<>A) then bugmes();quit;end
  460. v=sparse([],[],[6,1]);
  461. a1=a;a1(:,1)=v;A=full(a);V=full(v);A(:,1)=V;
  462. if or(full(a1)<>A) then bugmes();quit;end
  463. a1=a;a1(:,2)=v;A=full(a);V=full(v);A(:,2)=V;
  464. if or(full(a1)<>A) then bugmes();quit;end
  465. a1=a;a1(:,1)=v;A=full(a);V=full(v);A(:,1)=V;
  466. if or(full(a1)<>A) then bugmes();quit;end
  467. a1=a;a1(:,2)=v;A=full(a);V=full(v);A(:,2)=V;
  468. if or(full(a1)<>A) then bugmes();quit;end
  469. a1=a;a1(:,3)=v;A=full(a);V=full(v);A(:,3)=V;
  470. if or(full(a1)<>A) then bugmes();quit;end
  471. //
  472. a1=a;a1(:,7)=v;A=full(a);V=full(v);A(:,7)=V;
  473. if or(full(a1)<>A) then bugmes();quit;end
  474. a1=a;a1(:,8)=v;A=full(a);V=full(v);A(:,8)=V;
  475. if or(full(a1)<>A) then bugmes();quit;end
  476. b=sparse([],[],[6,2]);
  477. a1=a;a1(:,[1 3])=b;A=full(a);B=full(b);A(:,[1 3])=B;
  478. if or(full(a1)<>A) then bugmes();quit;end
  479. a1=a;a1(:,[3 1])=b;A=full(a);B=full(b);A(:,[3 1])=B;
  480. if or(full(a1)<>A) then bugmes();quit;end
  481. b=sparse([],[],[6,3]);
  482. a1=a;a1(:,[1 3 1])=b;A=full(a);B=full(b);A(:,[1 3 1])=B;
  483. if or(full(a1)<>A) then bugmes();quit;end
  484. a=sparse([],[],[6,6]);
  485. a1=a;a1(:,1)=v;A=full(a);V=full(v);A(:,1)=V;
  486. if or(full(a1)<>A) then bugmes();quit;end
  487. a1=a;a1(:,2)=v;A=full(a);V=full(v);A(:,2)=V;
  488. if or(full(a1)<>A) then bugmes();quit;end
  489. a1=a;a1(:,1)=v;A=full(a);V=full(v);A(:,1)=V;
  490. if or(full(a1)<>A) then bugmes();quit;end
  491. a1=a;a1(:,2)=v;A=full(a);V=full(v);A(:,2)=V;
  492. if or(full(a1)<>A) then bugmes();quit;end
  493. a1=a;a1(:,3)=v;A=full(a);V=full(v);A(:,3)=V;
  494. if or(full(a1)<>A) then bugmes();quit;end
  495. //
  496. a1=a;a1(:,7)=v;A=full(a);V=full(v);A(:,7)=V;
  497. if or(full(a1)<>A) then bugmes();quit;end
  498. a1=a;a1(:,8)=v;A=full(a);V=full(v);A(:,8)=V;
  499. if or(full(a1)<>A) then bugmes();quit;end
  500. b=sparse([1 2;2 1;3 1;3 2;4 1;6 2],10:15,[6,2]);
  501. a1=a;a1(:,[1 3])=b;A=full(a);B=full(b);A(:,[1 3])=B;
  502. if or(full(a1)<>A) then bugmes();quit;end
  503. a1=a;a1(:,[3 1])=b;A=full(a);B=full(b);A(:,[3 1])=B;
  504. if or(full(a1)<>A) then bugmes();quit;end
  505. b=sparse([1 2;2 1;2 3;3 1;3 2;4 1;5 3;6 2],10:17,[6,3]);
  506. a1=a;a1(:,[1 3 1])=b;A=full(a);B=full(b);A(:,[1 3 1])=B;
  507. if or(full(a1)<>A) then bugmes();quit;end
  508. v=sparse([],[],[6,1]);
  509. a1=a;a1(:,1)=v;A=full(a);V=full(v);A(:,1)=V;
  510. if or(full(a1)<>A) then bugmes();quit;end
  511. a1=a;a1(:,2)=v;A=full(a);V=full(v);A(:,2)=V;
  512. if or(full(a1)<>A) then bugmes();quit;end
  513. a1=a;a1(:,1)=v;A=full(a);V=full(v);A(:,1)=V;
  514. if or(full(a1)<>A) then bugmes();quit;end
  515. a1=a;a1(:,2)=v;A=full(a);V=full(v);A(:,2)=V;
  516. if or(full(a1)<>A) then bugmes();quit;end
  517. a1=a;a1(:,3)=v;A=full(a);V=full(v);A(:,3)=V;
  518. if or(full(a1)<>A) then bugmes();quit;end
  519. //
  520. a1=a;a1(:,7)=v;A=full(a);V=full(v);A(:,7)=V;
  521. if or(full(a1)<>A) then bugmes();quit;end
  522. a1=a;a1(:,8)=v;A=full(a);V=full(v);A(:,8)=V;
  523. if or(full(a1)<>A) then bugmes();quit;end
  524. b=sparse([1 2;2 1;3 1;3 2;4 1;6 2],10:15,[6,2]);
  525. a1=a;a1(:,[1 3])=b;A=full(a);B=full(b);A(:,[1 3])=B;
  526. if or(full(a1)<>A) then bugmes();quit;end
  527. a1=a;a1(:,[3 1])=b;A=full(a);B=full(b);A(:,[3 1])=B;
  528. if or(full(a1)<>A) then bugmes();quit;end
  529. b=sparse([1 2;2 1;2 3;3 1;3 2;4 1;5 3;6 2],10:17,[6,3]);
  530. a1=a;a1(:,[1 3 1])=b;A=full(a);B=full(b);A(:,[1 3 1])=B;
  531. if or(full(a1)<>A) then bugmes();quit;end
  532. v=sparse([],[],[6,1]);
  533. a1=a;a1(:,1)=v;A=full(a);V=full(v);A(:,1)=V;
  534. if or(full(a1)<>A) then bugmes();quit;end
  535. a1=a;a1(:,2)=v;A=full(a);V=full(v);A(:,2)=V;
  536. if or(full(a1)<>A) then bugmes();quit;end
  537. a1=a;a1(:,1)=v;A=full(a);V=full(v);A(:,1)=V;
  538. if or(full(a1)<>A) then bugmes();quit;end
  539. a1=a;a1(:,2)=v;A=full(a);V=full(v);A(:,2)=V;
  540. if or(full(a1)<>A) then bugmes();quit;end
  541. a1=a;a1(:,3)=v;A=full(a);V=full(v);A(:,3)=V;
  542. if or(full(a1)<>A) then bugmes();quit;end
  543. //
  544. a1=a;a1(:,7)=v;A=full(a);V=full(v);A(:,7)=V;
  545. if or(full(a1)<>A) then bugmes();quit;end
  546. a1=a;a1(:,8)=v;A=full(a);V=full(v);A(:,8)=V;
  547. if or(full(a1)<>A) then bugmes();quit;end
  548. b=sparse([],[],[6,2]);
  549. a1=a;a1(:,[1 3])=b;A=full(a);B=full(b);A(:,[1 3])=B;
  550. if or(full(a1)<>A) then bugmes();quit;end
  551. a1=a;a1(:,[3 1])=b;A=full(a);B=full(b);A(:,[3 1])=B;
  552. if or(full(a1)<>A) then bugmes();quit;end
  553. b=sparse([],[],[6,3]);
  554. a1=a;a1(:,[1 3 1])=b;A=full(a);B=full(b);A(:,[1 3 1])=B;
  555. if or(full(a1)<>A) then bugmes();quit;end
  556. // row column insertion
  557. a=sparse([1 6;1 5;1 3;2 4;2 1;4 4;4 3;5 1;6 6],1:9,[6 6]);
  558. b=sparse([1 2;1 3;3 3],-(1:3),[3,3]);
  559. a1=a;a1(1,1)=sparse(30);A=full(a);A(1,1)=30;
  560. if or(full(a1)<>A) then bugmes();quit;end
  561. a1=a;a1(1,6)=sparse(30);A=full(a);A(1,6)=30;
  562. if or(full(a1)<>A) then bugmes();quit;end
  563. a1=a;a1(1,8)=sparse(30);A=full(a);A(1,8)=30;
  564. if or(full(a1)<>A) then bugmes();quit;end
  565. a1=a;a1(1:3,1:3)=b;A=full(a);A(1:3,1:3)=full(b);
  566. if or(full(a1)<>A) then bugmes();quit;end
  567. a1=a;a1(1:3,6:8)=b;A=full(a);A(1:3,6:8)=full(b);
  568. if or(full(a1)<>A) then bugmes();quit;end
  569. a1=a;a1(6:8,1:3)=b;A=full(a);A(6:8,1:3)=full(b);
  570. if or(full(a1)<>A) then bugmes();quit;end
  571. a1=a;a1([3 2 1],1:3)=b;A=full(a);A([3 2 1],1:3)=full(b);
  572. if or(full(a1)<>A) then bugmes();quit;end
  573. a1=a;a1([1 2 1],1:3)=b;A=full(a);A([1 2 1],1:3)=full(b);
  574. if or(full(a1)<>A) then bugmes();quit;end
  575. a1=a;a1([3 2 1],[3 2 1])=b;A=full(a);A([3 2 1],[3 2 1])=full(b);
  576. if or(full(a1)<>A) then bugmes();quit;end
  577. a1=a;a1([1 2 1],[3 2 1])=b;A=full(a);A([1 2 1],[3 2 1])=full(b);
  578. if or(full(a1)<>A) then bugmes();quit;end
  579. a1=a;a1([1 2 1],[1 2 1])=b;A=full(a);A([1 2 1],[1 2 1])=full(b);
  580. if or(full(a1)<>A) then bugmes();quit;end
  581. //sparse full
  582. a=sparse([1 6;1 5;1 3;2 4;2 1;4 4;4 3;5 1;6 6],1:9,[6 6]);
  583. vt=11:16;
  584. // full line insertion
  585. //----------------------
  586. a1=a;a1(1,:)=vt;A=full(a);A(1,:)=vt;
  587. if or(full(a1)<>A) then bugmes();quit;end
  588. a1=a;a1(3,:)=vt;A=full(a);A(3,:)=vt;
  589. if or(full(a1)<>A) then bugmes();quit;end
  590. a1=a;a1(1,:)=vt;A=full(a);A(1,:)=vt;
  591. if or(full(a1)<>A) then bugmes();quit;end
  592. a1=a;a1(3,:)=vt;A=full(a);A(3,:)=vt;
  593. if or(full(a1)<>A) then bugmes();quit;end
  594. a1=a;a1(7,:)=vt;A=full(a);A(7,:)=vt;
  595. if or(full(a1)<>A) then bugmes();quit;end
  596. a1=a;a1(8,:)=vt;A=full(a);A(8,:)=vt;
  597. if or(full(a1)<>A) then bugmes();quit;end
  598. a1=a;a1([1 3],:)=[vt;2*vt];A=full(a);A([1 3],:)=[vt;2*vt];
  599. if or(full(a1)<>A) then bugmes();quit;end
  600. a1=a;a1([3 1],:)=[vt;2*vt];A=full(a);A([3 1],:)=[vt;2*vt];
  601. if or(full(a1)<>A) then bugmes();quit;end
  602. a1=a;a1([1 3 1],:)=[vt;2*vt;3*vt];A=full(a);A([1 3 1],:)=[vt;2*vt;3*vt];
  603. if or(full(a1)<>A) then bugmes();quit;end
  604. a=sparse([],[],[6,6]);
  605. a1=a;a1(1,:)=vt;A=full(a);A(1,:)=vt;
  606. if or(full(a1)<>A) then bugmes();quit;end
  607. a1=a;a1(3,:)=vt;A=full(a);A(3,:)=vt;
  608. if or(full(a1)<>A) then bugmes();quit;end
  609. a1=a;a1(1,:)=vt;A=full(a);A(1,:)=vt;
  610. if or(full(a1)<>A) then bugmes();quit;end
  611. a1=a;a1(3,:)=vt;A=full(a);A(3,:)=vt;
  612. if or(full(a1)<>A) then bugmes();quit;end
  613. a1=a;a1(7,:)=vt;A=full(a);A(7,:)=vt;
  614. if or(full(a1)<>A) then bugmes();quit;end
  615. a1=a;a1(8,:)=vt;A=full(a);A(8,:)=vt;
  616. if or(full(a1)<>A) then bugmes();quit;end
  617. b=[1:6;11:16];
  618. a1=a;a1([1 3],:)=b;A=full(a);A([1 3],:)=b;
  619. if or(full(a1)<>A) then bugmes();quit;end
  620. a1=a;a1([3 1],:)=b;A=full(a);A([3 1],:)=b;
  621. if or(full(a1)<>A) then bugmes();quit;end
  622. b=[1:6;11:16;21:26];
  623. a1=a;a1([1 3 1],:)=b;A=full(a);A([1 3 1],:)=b;
  624. if or(full(a1)<>A) then bugmes();quit;end
  625. // full column insertion
  626. //----------------------
  627. a=sparse([1 6;1 5;1 3;2 4;2 1;4 4;4 3;5 1;6 6],1:9,[6 6]);
  628. v=(1:6)';
  629. a1=a;a1(:,1)=v;A=full(a);A(:,1)=v;
  630. if or(full(a1)<>A) then bugmes();quit;end
  631. a1=a;a1(:,2)=v;A=full(a);A(:,2)=v;
  632. if or(full(a1)<>A) then bugmes();quit;end
  633. a1=a;a1(:,1)=v;A=full(a);A(:,1)=v;
  634. if or(full(a1)<>A) then bugmes();quit;end
  635. a1=a;a1(:,2)=v;A=full(a);A(:,2)=v;
  636. if or(full(a1)<>A) then bugmes();quit;end
  637. a1=a;a1(:,3)=v;A=full(a);A(:,3)=v;
  638. if or(full(a1)<>A) then bugmes();quit;end
  639. //
  640. a1=a;a1(:,7)=v;A=full(a);A(:,7)=v;
  641. if or(full(a1)<>A) then bugmes();quit;end
  642. a1=a;a1(:,8)=v;A=full(a);A(:,8)=v;
  643. if or(full(a1)<>A) then bugmes();quit;end
  644. b=[(1:6)' (11:16)'];
  645. a1=a;a1(:,[1 3])=b;A=full(a);A(:,[1 3])=b;
  646. if or(full(a1)<>A) then bugmes();quit;end
  647. a1=a;a1(:,[3 1])=b;A=full(a);A(:,[3 1])=b;
  648. if or(full(a1)<>A) then bugmes();quit;end
  649. b=[(1:6)' (11:16)' (21:26)'];
  650. a1=a;a1(:,[1 3 1])=b;A=full(a);A(:,[1 3 1])=b;
  651. if or(full(a1)<>A) then bugmes();quit;end
  652. //********
  653. // row column insertion
  654. a=sparse([1 6;1 5;1 3;2 4;2 1;4 4;4 3;5 1;6 6],1:9,[6 6]);
  655. b=-[1 2 3;4 5 6;7 8 9];
  656. a1=a;a1(1,1)=30;A=full(a);A(1,1)=30;
  657. if or(full(a1)<>A) then bugmes();quit;end
  658. a1=a;a1(1,6)=30;A=full(a);A(1,6)=30;
  659. if or(full(a1)<>A) then bugmes();quit;end
  660. a1=a;a1(1,8)=30;A=full(a);A(1,8)=30;
  661. if or(full(a1)<>A) then bugmes();quit;end
  662. a1=a;a1(1:3,1:3)=b;A=full(a);A(1:3,1:3)=b;
  663. if or(full(a1)<>A) then bugmes();quit;end
  664. a1=a;a1(1:3,6:8)=b;A=full(a);A(1:3,6:8)=b;
  665. if or(full(a1)<>A) then bugmes();quit;end
  666. a1=a;a1(6:8,1:3)=b;A=full(a);A(6:8,1:3)=b;
  667. if or(full(a1)<>A) then bugmes();quit;end
  668. a1=a;a1([3 2 1],1:3)=b;A=full(a);A([3 2 1],1:3)=b;
  669. if or(full(a1)<>A) then bugmes();quit;end
  670. a1=a;a1([1 2 1],1:3)=b;A=full(a);A([1 2 1],1:3)=b;
  671. if or(full(a1)<>A) then bugmes();quit;end
  672. a1=a;a1([3 2 1],[3 2 1])=b;A=full(a);A([3 2 1],[3 2 1])=b;
  673. if or(full(a1)<>A) then bugmes();quit;end
  674. a1=a;a1([1 2 1],[3 2 1])=b;A=full(a);A([1 2 1],[3 2 1])=b;
  675. if or(full(a1)<>A) then bugmes();quit;end
  676. a1=a;a1([1 2 1],[1 2 1])=b;A=full(a);A([1 2 1],[1 2 1])=b;
  677. if or(full(a1)<>A) then bugmes();quit;end
  678. // vector insertion
  679. v=sparse([1 1;3 1;6 1],[10 11 12],[6 1]);
  680. v1=v;v1(1)=33;V=full(v);V(1)=33;
  681. if or(full(v1)<>V) then bugmes();quit;end
  682. v1=v;v1(2)=33;V=full(v);V(2)=33;
  683. if or(full(v1)<>V) then bugmes();quit;end
  684. v1=v;v1(8)=33;V=full(v);V(8)=33;
  685. if or(full(v1)<>V) then bugmes();quit;end
  686. v1=v;v1([1 2 8])=[5;10;33];V=full(v);V([1 2 8])=[5;10;33];
  687. if or(full(v1)<>V) then bugmes();quit;end
  688. v1=v;v1(:)=[];
  689. if or(full(v1)<>[]) then bugmes();quit;end
  690. v1=v;v1(1)=sparse(33);V=full(v);V(1)=33;
  691. if or(full(v1)<>V) then bugmes();quit;end
  692. v1=v;v1(2)=sparse(33);V=full(v);V(2)=33;
  693. if or(full(v1)<>V) then bugmes();quit;end
  694. v1=v;v1(8)=sparse(33);V=full(v);V(8)=33;
  695. if or(full(v1)<>V) then bugmes();quit;end
  696. v1=v;v1([1 2 8])=sparse([5;10;33]);V=full(v);V([1 2 8])=[5;10;33];
  697. if or(full(v1)<>V) then bugmes();quit;end
  698. v1=v;v1([1 2 1])=sparse([5;10;33]);V=full(v);V([1 2 1])=[5;10;33];
  699. if or(full(v1)<>V) then bugmes();quit;end
  700. v1=v;v1(:)=[];
  701. if or(full(v1)<>[]) then bugmes();quit;end
  702. v1=v;v1(:)=sparse([2 1],44,[6 1]);V=full(v);V(:)=[0;44;0;0;0;0];
  703. if or(full(v1)<>V) then bugmes();quit;end
  704. v=v';
  705. v1=v;v1(1)=33;V=full(v);V(1)=33;
  706. if or(full(v1)<>V) then bugmes();quit;end
  707. v1=v;v1(2)=33;V=full(v);V(2)=33;
  708. if or(full(v1)<>V) then bugmes();quit;end
  709. v1=v;v1(8)=33;V=full(v);V(8)=33;
  710. if or(full(v1)<>V) then bugmes();quit;end
  711. v1=v;v1([1 2 8])=[5 10 33];V=full(v);V([1 2 8])=[5 10 33];
  712. if or(full(v1)<>V) then bugmes();quit;end
  713. v1=v;v1(1)=sparse(33);V=full(v);V(1)=33;
  714. if or(full(v1)<>V) then bugmes();quit;end
  715. v1=v;v1(2)=sparse(33);V=full(v);V(2)=33;
  716. if or(full(v1)<>V) then bugmes();quit;end
  717. v1=v;v1(8)=sparse(33);V=full(v);V(8)=33;
  718. if or(full(v1)<>V) then bugmes();quit;end
  719. v1=v;v1([1 2 8])=sparse([5 10 33]);V=full(v);V([1 2 8])=[5 10 33];
  720. if or(full(v1)<>V) then bugmes();quit;end
  721. v1=v;v1([1 2 1])=sparse([5 10 33]);V=full(v);V([1 2 1])=[5 10 33];
  722. if or(full(v1)<>V) then bugmes();quit;end
  723. v1=v;v1(:)=sparse([1 2],44,[1,6]);V=full(v);V(:)=[0 44 0 0 0 0];
  724. if or(full(v1)<>V) then bugmes();quit;end
  725. v1=v;v1(1)=[];V=full(v);V(1)=[];
  726. if or(full(v1)<>V) then bugmes();quit;end
  727. //test des comparaisons
  728. a=sparse([1 6;1 5;1 3;2 4;2 1;4 4;4 3;5 1;6 6],1:9,[6 6]);
  729. b=sparse([1 6;1 2;6 5],[10;-1;-1],[6 6]);
  730. if full(a==a)<>full(a)==full(a) then bugmes();quit;end
  731. if full(a<>a)<>(full(a)<>full(a)) then bugmes();quit;end
  732. if full(a>sparse(5))<>(full(a)>5) then bugmes();quit;end
  733. if full(sparse(5)>a)<>(5>full(a)) then bugmes();quit;end
  734. if full(b>a)<>(full(b)>full(a))  then bugmes();quit;end
  735. if full(a==full(a))<>full(a)==full(a) then bugmes();quit;end
  736. if full(a<>full(a))<>(full(a)<>full(a)) then bugmes();quit;end
  737. if full(a>5)<>(full(a)>5) then bugmes();quit;end
  738. if full(5>a)<>(5>full(a)) then bugmes();quit;end
  739. if full(b>full(a))<>(full(b)>full(a))  then bugmes();quit;end
  740. if full(full(a)==a)<>full(a)==full(a) then bugmes();quit;end
  741. if full(full(a)<>a)<>(full(a)<>full(a)) then bugmes();quit;end
  742. if full(full(a)>sparse(5))<>(full(a)>5) then bugmes();quit;end
  743. if full(full(b)>a)<>(full(b)>full(a))  then bugmes();quit;end
  744. a=sparse([1 1;3 1;6 1],[10 11 12],[6 1]);
  745. if full(a==a)<>full(a)==full(a) then bugmes();quit;end
  746. if full(a<>a)<>(full(a)<>full(a)) then bugmes();quit;end
  747. if full(a>sparse(5))<>(full(a)>5) then bugmes();quit;end
  748. if full(sparse(5)>a)<>(5>full(a)) then bugmes();quit;end
  749. if full(a==full(a))<>full(a)==full(a) then bugmes();quit;end
  750. if full(a<>full(a))<>(full(a)<>full(a)) then bugmes();quit;end
  751. if full(a>5)<>(full(a)>5) then bugmes();quit;end
  752. if full(5>a)<>(5>full(a)) then bugmes();quit;end
  753. if full(full(a)==a)<>full(a)==full(a) then bugmes();quit;end
  754. if full(full(a)<>a)<>(full(a)<>full(a)) then bugmes();quit;end
  755. if full(full(a)>sparse(5))<>(full(a)>5) then bugmes();quit;end
  756. a=a';
  757. if full(a==a)<>full(a)==full(a) then bugmes();quit;end
  758. if full(a<>a)<>(full(a)<>full(a)) then bugmes();quit;end
  759. if full(a>sparse(5))<>(full(a)>5) then bugmes();quit;end
  760. if full(sparse(5)>a)<>(5>full(a)) then bugmes();quit;end
  761. if full(a==full(a))<>full(a)==full(a) then bugmes();quit;end
  762. if full(a<>full(a))<>(full(a)<>full(a)) then bugmes();quit;end
  763. if full(a>5)<>(full(a)>5) then bugmes();quit;end
  764. if full(5>a)<>(5>full(a)) then bugmes();quit;end
  765. if full(full(a)==a)<>full(a)==full(a) then bugmes();quit;end
  766. if full(full(a)<>a)<>(full(a)<>full(a)) then bugmes();quit;end
  767. if full(full(a)>sparse(5))<>(full(a)>5) then bugmes();quit;end
  768. a=sparse([1 6;1 5;1 3;2 4;2 1;4 4;4 3;5 1;6 6],1:9,[6 6]);
  769. b=sparse([1 6;1 2;5 5],[10;-1;-1],[5 6]);
  770. if a==b<>%f then bugmes();quit;end
  771. if a<>b<>%t then bugmes();quit;end
  772.